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NOMENCLATURE 



2 / 

0 thermal diffusivity, cm /sec 

a coefficient of linear thermal expansion, l/°K 

A relative amplitude of electromagnetic radiation 

b y-coordinate of vertex- of parabola 

c x-coordinate of vertex of parabola 

d thickness of metal sheet or mirror, cm 

D diameter of mirror, meters 

f ratio of aperture to focal length 

F focal length, meters 

Fo Fourier number, nondimensional time 

2 2 

1 intensity, watts/cm or kW/cm 

2 2 

intensity of laser beam, watt s/cra , kW/cm 
k thermal conductivity, watts/cm °K 

I length of object expanding 

n unit vector normal to surface 

p distance from focus to directrix, meters 

r reflectivity, or radial coordinate on mirror, meters 

R characteristic dimension of laser beam, or nondimensional mirror radius 

t temperature, °K 

t^ characteristic temperature, % 

t^ initial temperature, 

T nondimensional temperature 
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0 time, seconds; or angle between beam aoid normal, degrees 

X coordinate for distance through sheet nf metal or mirror, centimeters, 

meters 

1 nondiroensional distance through plate, x/d 
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I. INTRODUCTION 



Among the applications of high energy lasers is welding. For this applica- 
tion and others, one would like an estimate of the time required to achieve a 
specified temperature. Other information that is useful is temperature 
distribution within the metal. 

A set of graphs has been developed which permits rapid evaluation of time- 
to -temperature, temperature profiles, etc. First, the theory is given for the 
problem. Next, the limitations are outlined clearly. Finally, the graphs are 
discussed and sample problems are woiiced. 



Equation (l) is solved for a plate of thickness d with a constant heat flxtx at 
the surface x“ 0, i. e., the front wall. The back wall, x » d, is assumed to 
be insulated so that 



II. THEORY AND ANALYSIS 



The time-dependent, one -dimensional equation for heat conduction is 



9t ^ ^ 

86 8x^ 



( 1 ) 




( 2 ) 



2 

The constant heat flux is due to a laser beam of intensity I, watts/cm . 
Consequently, the boundary condition on the front wall is 




( 3 ) 
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( 4 ) 



Introduce the following nondi»enBional Tariables into equation* (l) and (3) 
(a), I ■ x/d} (b), T = (t - t^Vt^; (c), t^ » IdA? (d) Fo »a ©/d^ 
giTing 



8T ^ sfr 

»Fo 8X^ 



(5) 



and 




( 6 ) 



The solution to the problea formulated above ie 



T - Fo ♦ 



1 2 -FoTT^n^ 

3 2 ^ n 



cos mX 



(7) 



n*l 

Due to the exponential term in the series, for Fo > 0,4 the error in dropping 
the series is less than 1 per cent. 



III. LIMITATIONS OF ANALYTICAL M3DEL 



The assumption of one-dimensional heating is valid so long as R > > d, \diere 
R is a characteristic dimension of the laser beam. The assumption of constant 
flux into the metal, as stated by equation (3), aay or may not be realistic. 

This point will be discussed in more detail below. 

Implicit in the equations (l) to (7) is the fact that material properties 
remain constant. Thermal conductivity, k, is a function of temperature. Like- 
wise thermal diffusivity, a , is a function of temperature. 

Ret\im again to the question of constant heat flttx due to laser radiation. 
The laser beam has the intensity I^, The value of intensity that is correct for 
equation (3) is 

I » I]j®os 0(1 - r) (8) 
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Equation (8) accounts for the spreading of the radiation due to angle 0 
(See Fig. 1. ) emd the reflectivity. The reflectivity, r, is a weak function 
of temperature. 

The charts are of value only up to the melting temperature. All 
parameters in the problem vary rapidly with temperature as phase changes from 
solid to liquid. Furthermore, the heat of vaporiaation is not incorporated in 
the model. When the surface becomes molten, vapor forms. The vapor may be 
absorbing} as a consequence, the intensity of radiation may be greatly reduced. 

In suimnary, the model provides adequate answers up to the melting tempera- 
ture if the radiant heat flux, I, is correctly deteimiined by equation (7), 

IV. DISCUSSION OF GRAPHS 

Figures 2 to 4 are graphs which have been calculated using equation (7) 

and the appropriate nondimens ional variables, equations (4a) to (4d), which 

were defined earlier. Note that there are three graphs, i. e. , one in 

upper left hand (ULH) comer, one in upper right hand (URH) comer, and one in 

lower right hand (LRH) corner. The primary graph is in the URH corner; both 

ULH and LRH graphs are for convenience in the use of the primary graph. The 

URH graph is a plot of nondimens ional temperature, T, in terms of non- 

dimensional time, Fo. The Fourier number is essentially a nondimens ional 

time. To use the graphs, one must know the intensity, I; the thickness of the 

plate, d; the thermal conductivity, k; and thermal diffusivity, a . A 

characteristic temperature, t , is calculated using equation (4e); t is used 

c c 

to enter the ULH graph. 

The LRH graph has been drawn for aluminum. It is a plot of 
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Figure 1. Geometry of Heating by Laser Radiation. 
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2 

For alumlmun» the thermal diffusivity is 0.521 cm /sec. The time, 6, Is used 
as a parameter in the LRH graph. 

Both the ULH and URH graphs eo'e universal; i.e., they can be applied to any 
homogeneous substance. However, the LRH graph is applicable only for aluminum. 

A simple calculation can be used to convert the time, 0, read from the graph 
to the time, 0', for other materials. This calculation is based on equation (9) 
as follows: 

0' * 0a /a ' 

The other material has thermal diffusivlty, a'. For example, suppose 
2 

* 0,2605 cm /sec for material "I"; hence 

e« - (1»2 second)(0.521.ca^/sec} ^ 3.4 second 
(0.2605 cm /sec) 

With the simple calculation outlined above, the graphs of Figs. 2, 3, and 4 
are universal. 

Some sample problems will now be discussed. 

A. TEMPERATUEE AT A GIVEN TIME AND LOCATION 
Sample Problem 1 

Assume a plate thickness of d « 1.3 cm and an initial temperature of 

tp * 300°K. Assume I = 1000 watt/cm^ and k *= 1.3 watt/cm°K. If time 

0 » 0.4 second has elapsed since the laser beam was switched on, what is the 

temperature at front face? 

Calculate t . 

c 

= 1000°K 



. (1000)(1.3) 
^c (1.3) 



Enter the LRH graph with d » 1.3 ca. Move to the right to the curve for 

6 * 0.4. Then move upward to the curve for X » 0, which is the front face. 

Note that Fo = 0.12 and T = .39. Jfove toward the left until the line for 

t = 1000*^K is intersected. Then move down to the axis and read t - t~ = 390, 

c 0 

Since the initial temperature was 300*^K, the front face has a temperature 
t = 690°K. 

Sample Prohlen 2 

Sample Problem 2 is the same as Problem 1 except that the temperature at 
X ® 0.52 cm is desired. Calculate X = x/d = 0,52/1.3 * 0.4. Where the 
vertical lino running between the UIH and URH graphs intersects X = 0.4 curve, 
draw a horizontal line. At the intersection with the t = 1000°K curve, 
descend vertically to the axis. Road t - t^ » 115* At a distance x = 0.52 ca 
and a time 0 » 0.4 sec, the temperature is 415°K. 

B. TIME TO MELT FRONT SURFACE 
Sample Problem 3 

The melting temperature of aluminum is 933°K. Assume t^ * 333°K. 

2 

If a beam has an Intensity of 2000 watt/cm , how much time is required to 

begin melting the front surface of an alximlnum plate with thickness 0.65 cm? 

Calculate t . 

c 

. 4000“k 

Subtracting 333°K from 933°K, one enters the ULH graph at 600°K and moves 

upward to t * 4000°K curve. Moving horizontally, one crosses the T-axis 
c 

at 0,152 and stops at the intersection of X = 0 curve. Then a line is dropped 
vertically crossing the Fourier number axis at 0.018. Draw a horizontal curve 
in the LRH graph at d ■ 0,65 cm. The Intersection of the horizontal and vertical 
lines shows that 0 = .015 second. Note that the front face is ready to melt while 
the tenqjerature of the back face (X * l) remains t^ ■ 333°K. 
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V. MIBROR DISTORTION DDE TO RADIANT HEATING 



High powered lasers often use metal mirrors which undergo thermal expansion 
a«i optical distortion due to absorption of radiation. Depending on mirror 
geometry, complex, three-dimensional calculations may be required. In this section, 
a sample calculation for determining one-dimensional mirror distortion is Illustrated, 
The formula for thermal one-dimensional strain is 

^ « AAt (10) 

where a is the coefficient of linear thermal expansion, t is the length of 

the object expanding, A-t is the increase in its length, and At is the temperature 

change involved. 

The linear thermal expansion coefficient, a, is a function of temperature. 

Table I gives its values for different temperatures of almainum. Values for 
other metals as a function of temperature can be found in the American Institute 
of Physics Handbook ^ef. 



Table I. Temperature Variation of Linear Thermal 

Expansion Coefficient for Aluminum /Ref. 1^. 



t (°K) 


a (10"^/°K) 


300 


23.2 


350 


24.1 


400 


24.9 


500 


26.4 


600 


28.3 


700 


30.7 
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The expansion is given hy equation (11) 



d 



/ 



a(t)(t - tQ)dx 



( 11 ) 



where x is the coordinate for the distance through the mirror and d is the 
mirror thickness. For this sample problem, numerical evaluation has been 
chosen. The following characteristics for the mirror were used: parabolic 

2 

mirror, mirror radius * 0.05 m, f number = 1.0, intensity of beam = 40 kW cm , 
mirror absorption » 0,01, and Fresnel number * 10, 

Using the facts that the mirror is parabolic with a f number of 1.0 and 
thickness at the center of 2 cm, the thickness of the mirror at other radial 
points can be calculated, 

^ Focal length ^ F 
Diameter D 

F * Uf = (0.1 m)(l) = 0.1 m 
The standard form for a parabola is 

2p(y - b) * (r - c)^ (12) 

where p is the distance from the directrix to the focus of the parabola and 
(c,b) are the coordinates of its vertex. Thus p/2 is the focal length of the 
mirror. So p/2 ® 0,1m, and p = 0,2 m. Putting the vertex (c,b) of the 
parabola at the origin yields c * b » 0, Therefore, the equation of the 
parabola is 2(0,2)(y - O) * (r - 0) or y * 2. 5r , Solving this equation for 
different values of r gives the increase in the thickness of the mirror at 
different radii. The actual thickness can be found ly adding 2 cm, which is the 
thickness at the center of the mirror. The thickness of the mirror at different 
radial locations is given in Table II. 
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Table II. Mirror Thickness, Intensity of Absorbed Power, 
and Characteristic Temperature 



R 


r, m 


d, cm 


2 

Intensity, w/cm 


t , °K 


0 


0.0 


2.000 


400 


615 


0.3 


0.015 


2.056 


296 


468 


0.7 


0.035 


2.306 


77 


137 


1.0 


0.05 


2.625 


5 


9.77 



To find the Intensity of the laser radiation at Tarious radii on the mirror, 

the amplitude profile from Figure 3.11 from Lengyel ^ef. ^ was used. The absorbed 

2 

intensity at the center of the mirror was assumed to be 400 watts/cm which 

was obtained by multiplying radiation intensity fcy mirror absorption 
2 

(40kW/cm X O.Ol). Amplitude and intensity are related by 



h 

I« 




(13) 



where Aq and A^ are the amplitudes at two different points (from the referenced 
graph), and Iq is the known intensity at the mirror uis. The intensities also 
are tabulated in Table II. 



In this problem, one*<dimen8lonal heating is assumed; l.e. , the effects of 

heat transfer in the radial direction are ignored. To solve the problem, it is 

necessary to calculate t^. This can be done using equation (4c) where 

h = 1.3 watt/cm ^K. The values of t for different values of d are listed in 

c 

Table II. 



19 



To evaluate the integral numerically, the mirror is imagined as being 
divided into slabs. The expansion of each of these slabs is calculated 
based on the mean slab temperature. These calculations can be made using 
the set of graphs. Figs, 2 to 4, explained earlier in this paper. The 
calculation methods become apparent by referring to Table III which svimitiarizes 
the results for the centerline position on the mirror. 

The graph. Fig. 5, is a profile of the thermal expansion. From the 
graph, it can be seen that considerable distortion will occur, even at times 
as small as 0.3 sec. Modest quality optics have a tolerance less than one 
quarter of the wavelength of the light. The distortion in this problem 
would be considerably more than a quarter wavelength. 

VI. CONCLUSIONS 



Graphs have been developed for the solution of many types of problems 
involved with laser heating. Rapid solutions become possible by suitable 
use of the figures. Sample problems were presented for temperature profiles 
within the plate, time to melt the front surface, and mirror distortion due to 
thermal expansion. 



Table III. 


Summar 7 of 


Calculations for Thermal Expansion at Mirror Center 


X = x/d 


X, cm 




<t - tp> a 


Ax 


8x, microns 


0 


0 


uo 














no .252 


.2 


5.544 


0.1 


0.2 


80 














65 .244 


.2 


3.172 


0.2 


0.4 


50 














37.5 .239 


.2 


1.793 


0.3 


0.6 


25 














17.5 .235 


.2 


0.823 


0.4 


0.8 


10 














6 .233 


.4 


0.559 


0.6 


1.2 


2 














1 .232 


.4 


0.093 


0.8 


1.6 


0 














0 


- 


- 


1.0 


2.0 


0 












Total 


Thermal Expansion, tX 




11.983 microns 


^Obtained from Fig. 3 


with d = 2 


cm and t = 615°K. 
c 







< t - 1 q> average slab temperature in excess of initial temperature t^, °K. 
a coefficient of linear thermal expansion, microns/cm °K, 

Ax slab thickness, cm. 

8x thermal expansion of slab, microns. 
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Figure 5. Profile of Thermal Expansion at 300 Millisecs. 
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